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Abstract 

It is shown how the QCD sum rules can be applied for the investigation 
of the density dependence of the nucleon parameters. These characteristics 
can be expressed through the expectation values of QCD operators in nuclear 
matter. In certain approximations the expectation values are related to the 
observables. First applications of the approach reproduced some of the basic 
features of nuclear physics, providing also a new knowledge. The program of 
the future work is presented. The difficulties of the approach are discussed. 
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4. Summary 

The lecture is not addressed to the experts. The aim of the talk is rather 
to attract attention of the researchers who just started to study the sub- 
ject. The third section is based on the results obtained in collaboration with 
M. G. Ryskin, V. A. Sadovnikova and E. M. Levin. 

1 Motivation 

The theory of nuclear matter leaves some room for the improvement. This 
concerns the low densities, i.e. those close to the saturation value as well as 
the higher densities. 

Since the pioneering paper of Walecka pQ the partially successful Schrodinger 
phenomenology was succeeded by still more successful Dirac phenomenology. 
A nucleon in nuclear matter is treated as moving in superposition of the scalar 
and vector fields which are of several hundreds MeV. In the meson-exchange 
picture of nucleon interactions these fields originate from the exchange by o 
and uj mesons. This model is known as quantum hadrodynamics (QHD). It 
is quite successful in describing most of the properties of nucleons in both 
nuclear matter and finite nuclei [21 E] . However the model is not fundamen- 
tally complete. The weak points of QHD were reviewed by Negele j3] and 
by Sliv et al [S]. Here I mention that "cr-meson" is rather an effective way of 
describing the scalar interactions. Also the masses of a- and cu-mesons 
are so large that the scalar and vector interactions take place at the dis- 
tances where the nucleons cannot be treated as the point particles any more. 
Finally, the coupling constants g a and g u are the free parameters of QHD. 
They are chosen usually to fit the saturation properties. 

Thus it would be desirable to develop the approach which has the attrac- 
tive features of QHD, avoiding, however, the "meson-exchange" conception. 
This would enable us to avoid the controversy of the small distance descrip- 
tion mentioned above. 
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Another point is connected with the high densities. There are several 
interesting phenomena. The chiral phase transition and the disintegration of 
the nuclear matter to the quark-gluon plasma are much discussed nowadays. 
The possibility of the existence of the other phase states of nuclear matter 
containing the admixture of heavier baryons or of the "pion condensate" 
have been studied long ago [UJ [7] . These effects can be important for the 
astrophysics. However the QHD parameters g s ^ are defined at the saturation 
point. Hence, this approach cannot be expanded to the higher densities in a 
straightforward way. It is desirable to express the interactions in the scalar 
and vector channels through the observables whose density dependence can 
be found separately. This would enable to study the high density nuclear 
physics. 

There are chances that both requirements can be realized. Recall that the 
QCD sum rules (SR) method invented by Shifman et al. [Hj succeeded in ex- 
pressing the static properties of the hadrons through the vacuum expectation 
values of several simplest operators of the quark and gluon fields (QCD oper- 
ators of the lowest dimensions). If we succeed in expanding the SR approach 
for the case of finite densities, the in-medium modification of the values of 
the hadron parameters would be expressed through the in-medium values of 
the QCD condensates. Such approach would not require the conception of 
heavy meson (a and u) exchange. The calculation of the density dependence 
of QCD condensates would enable to use the approach in the broad interval 
of the density values. 

Now I give a brief review of the SR method in vacuum, focusing on the 
points which we shall need at the finite values of density. There are several 
detailed reviews of the SR approach P ITU]. 

2 QCD sum rules in vacuum 
2.1 Dispersion relations 

The basic point of SR method is the dispersion relation for the function G(q 2 ) 
which describes the propagation of the system with the quantum numbers of 
a hadron. In the simplest form it is 




(1) 
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In quantum mechanics G(q 2 ) is just the particle propagator. In the field 
theory different degrees of freedom are convenient in the different regions of 
the values of q 2 . In particular, for the system with the baryon and electric 
charges equal to unity Q — B — 1, the imaginary part ImG(fc 2 ) = at k 2 < 
m 2 with m being the position of the lowest lying pole, i.e. m is the proton 
mass. There are the other singularities at larger values of k 2 . There are the 
cuts corresponding to the systems "proton+pions" , etc. On the other hand, 
one can consider the system as that of three strongly interacting quarks. Such 
description becomes increasingly simple at q 2 — > — oo due to the asymptotic 
freedom of QCD. This means that at q 2 — > — oo the function G(q 2 ) can 
be presented as the power series of q~ 2 (and of QCD coupling constant a s ). 
The coefficients of the expansion are the expectation values of local operators 
constructed of quark and gluon fields. These expectation values are called 
"condensates". Thus such presentation known as operator power expansion 
(OPE) [TT], provides the perturbative expansion of the short-distance effects, 
while the nonperturbative physics is contained in the condensates. 

Technically this means that one should start with the general presentation 

G{q 2 ) = i J d 4 xe i{qx) (0\T{r](x)rtO)}\0) (2) 

with i] being the local operator with the quantum numbers of the considered 
system. The operator 77 (x) is the composition of the quark fields ip(x). For 
each quark field one can write in the lowest order of a s (forgetting for a while 
about colours) 

, im, 2 i 

(01^(^(0)10) = -V aP -- (0| : g(0)r^g(x) : [0) (3) 

with x = x^7 M . This is the direct consequence of Wick theorem. In Eq. 
(3) a and (3 are the Lorentz indices, m q stands for the quark mass. In the 
second term of the rhs of Eq. (3) T x is the complete set of the basic Dirac 
4x4 matrices with the scalar, vector, pseudoscalar, pseudovector and tensor 
structures. The first term in the rhs of Eq. (3) is just the free propagator. 

In the theories with the empty vacuum, e.g. in quantum electrodynamics 
the second term in rhs of Eq. (3) vanishes due to the normal ordering. In 
any field theory all the structures except the scalar one vanish due to the 
Lorentz invariance. In QCD the scalar term survives due to the spontaneous 
breaking of the chiral invariance. Thus 

(0|TgS(x)^(0)|0> = _L + y Z 2 * 2 got _ _L (0^(0)^)10) (4 ) 
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with "a" and "6" standing for the colour indices. Note that Eq.(4) has a 
simple physical meaning. The quark can propagate between the space-time 
points "0" and " X ctS cl free particle or by the exchange with the vacuum sea 
of the quark-antiquark pairs. Note, however, that Eqs. (3) and (4) are not 
presented in a gauge-invariant way (the operator q(x) depends on the gauge 
of the gluon fields). The rigorous and gauge- invariant form of Eq. (4) is 

<0|Tg£(s)#(0)|0> = _i_ jog + *K)/2* 2 ^ _ l_ m a {Q)q{m + 0(g2) 

(5) 

with the expansion in powers of x 2 corresponding to the expansion of G(q 2 ) 
in powers of q~ 2 . 

To demonstrate the power of the dispersion relations I present, following 
[TUj . the derivation of the well-known Gell-Mann, Oakes-Renner relation 
(GMOR) H21 

2 f 2 m 2 

(0\uu + dd\0) = - U n (6) 

TTl u -\- Vdfj. 

with f n and being the decay constant and the mass of the pion. Recall 
that Eq. (6) is true in the chiral limit m 2 —>■ 0, m q — > 0. 

The quantum numbers of pion can be carried by the axial current A^(x) = 
u(x) , j fl r y 5 d(x) as well as by the pseudoscalar current P(x) = iu(x)j 5 d(x) . 
Consider the dispersion relation for the function 

. q»fd*xe^(0\TA^x)P(0)\0) 
G{q) = i (7) 

using Eq. (5) for the quark propagators. The corresponding integral diverges 
at small x. Introducing a cutoff x 2 > L 2 one finds in the limit a s = 



8tt z —q z q z 



To obtain the rhs of Eq. (1) recall about the partial conservation by axial 
current (PCAC) expressed by the equation D^A^x) = ^/2f 7T m 2 r (p(x) with ip 
standing for the pion field. Present 

Im G(k 2 ) = (0|A M P|7r(A;))(7r(A;)|P|0)5(A; 2 -m^) + J R(A; 2 ) (9) 

with the term R(k 2 ) describing the higher lying states. Using PCAC one 
finds the term R(k 2 ) to contain one more factor m 2 compared to the first 
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term in rhs of Eq. (9). Hence, in the chiral limit we can neglect R(k 2 ) as 
well as the first in rhs of Eq. (8). Thus the dispersion relation is 



. (0\uu + dd\0) (0|A^|7r)(7r|P|0) , in , 

1 2 = 2 2 ' 10 

qz m z _ gz 

with k 2 = m 2 . Calculating the matrix elements in the rhs by using PCAC 
and assuming — q 2 ^> m 2 we come to Eq. (6). 

This derivation of GMOR relation is complementary to the standard one 
presented in the QCD textbooks (see, e.g. [13]) which is based on physics 
of small momenta q. Usually GMOR is treated as the way to determine the 
value of (0\uu + dd\0). However in the framework of the developed approach 
it can be viewed as the relation which expresses the combination of the pion 
parameters m 2 f 2 through the expectation value (0\uu + dd\0). 

2.2 Sum rules 

Unfortunately, the example considered above is the only case when the dis- 
persion relation takes such a simple form. Usually there is no reason to 
neglect the higher lying physical states with respect to the lowest one. If the 
second lowest singularity is the cut starting at the point W 2 h we can present 

Im G(k 2 ) = X 2 5(k 2 - m 2 ) + f(k 2 )6(k 2 -W 2 h ) (11) 

with A 2 being the residue while f(k 2 ) is the spectral function. Following 
previous discussion the lhs of Eq. (1) can be expanded in powers of q~ 2 and 
thus Eq. (1) takes the form 

G PE{q 2 ) = ^ + - / f{k2)dk2 (12) 
" m? — q 2 Ti J k 2 — q 2 

with the unknown parameters m, A 2 , and the unknown spectral function 
f(k 2 ). The aim of the SR is to obtain the parameters of the lowest lying 
state. Hence, the second term of rhs of Eq. (12) is treated approximately. 
The approximation is prompted by the asymptotic behaviour 

f{k 2 ) = ^AGoP E (k 2 ) (13) 
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at k 2 ^> \q 2 \. The integral over k 2 3> \q 2 \ provides the terms ~ ln(L 2 / — q 2 ) 
exceeding the contribution of the pole by this factor. The standard ansatz 
consists in extrapolation of Eq. (12) to the lower values of k 2 , replacing also 
the physical threshold W 2 h by the unknown effective threshold W 2 , i.e. 

w\ w 2 
and thus the dispersion relation (1) takes the form 

n ( 2\ ^ 2 1 f ^G PE{k 2 ) „ 2 / nK \ 

G ° Mq) = ^7 2 + 2^J k 2 -q 2 dk ■ (15) 

w 2 

Such approximation of the spectrum is known as the "pole + continuum" 
model. 

The lhs of Eq. (15) contains the QCD condensates. The rhs contains 
three unknown parameters m, X 2 and W 2 . However both lhs and rhs depend 
on q 2 . The OPE becomes increasingly true at large values of — q 2 . The "pole 
+ continuum" model has sense only if the contribution of the continuum, 
treated approximately does not exceed the contribution of the pole, treated 
exactly. Thus the model becomes increasingly true at small values of \q 2 \. 
The problem is to find the region of \q 2 \ where both OPE and "pole + 
continuum" model are valid. 

Such region is unlikely to exist in any channel of the dispersion relations 
presented by Eq. (15) with the necessary subtractions. To improve the 
overlap between the QCD and phenomenological descriptions Borel transform 
defined as 

B ^ - J^TrW'W*^ (16) 

Q 2 = -q 2 ; M 2 = Q 2 /n 

was used in jHj. There are several useful features of the Borel transform. 
It removes the divergent terms in the lhs of Eq. (15) which are caused by 
the free quark loops — see, e.g., Eq. (8). This happens, since the Borel 
transform eliminates all the polynomials in q 2 . Thus we do not need to make 
subtractions. Also it emphasizes the contribution of the lowest lying states 
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in rhs of Eq. (15) due to the relation 



B 



1 



,2/M 2 



(17) 



Q2 _j_ m 2 



= e 



The Borel transformed dispersion relations 



G OPE {M 2 ) = AV m2 / M2 + 



1 



/ 



dk 2 e~ k2 / M2 -AGoPEik 2 ) 



(18) 
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have been analyzed first for the vector mesons [B] and for the nucleons [Tlj . 
In both cases the regions of the values of the Borel mass M were found, where 
the matching of rhs and lhs of Eq. (18) was achieved. In particular, Ioffe 14 
calculated the value of the meson mass as the function of QCD condensates. 
Several terms of OPE appeared to be needed to provide the matching of rhs 
and lhs of Eq. (18). The result of ^3] can be treated as 



with C < being the function of the gluon condensate and of the four-quark 
condensate (0|gggg|0). The latter can be viewed as the expansion of the two- 
quark propagator, similar to Eq. (3). One can see that Eq. (19) has a simple 
physical meaning. The nucleon mass is caused by the quark exchange with 
the vacuum sea of qq pairs. The mechanism resembles that of the Nambu 
and Jona-Lasinio model. 

The QCD SR method was applied successfully to calculation of the static 
properties of mesons jHj and nucleons jT^l [T3J. It provided new knowledge 
as well. For example, the value of gluon condensate go = (0|— G^G^O) 
was extracted by Vainshtein et al. (THj from the analysis of leptonic decays 
of p and (p mesons and from QCD analysis of charmonium spectrum. This 
condensate is a very important characteristics of QCD vacuum, since it is 
directly related to the vacuum energy density. The investigations based on 
the vacuum QCD sum rules are going on until now. The latest HEP preprint 
17J which the value of go has been calculated more accurately was published 
several months ago. 



m = C (0\uu\0) 



(19) 
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3 QCD sum rules in nuclear matter 



3.1 The problems 

Now we discuss the possibility of the extension of QCD SR approach to the in- 
vestigation of the characteristics of the nucleons in nuclear matter. If we suc- 
ceed, the modification of the characteristics of the nucleon will be expressed 
through the expectation values of QCD operators in medium. Although some 
of the qualitative results may find the applications in the investigations of 
the finite nuclei, only the infinite nuclear matter will be considered below. 
In other words, the density of the distribution of the nucleons is the same in 
all the space points. 

Since the Lorentz invariance is lost, the correlation function in medium 



depends on two variables, but not on q 2 only. The spectrum of the function 
G m (q) is much more complicated, than that of the vacuum correlator G(q 2 ) 
defined by Eq. (2). The singularities can be connected with the nucleon 
(proton) placed into the matter, as well as with the matter itself. One of the 
problems is to find the proper variables, which would enable us to focus on 
the properties of our probe proton. 

In the papers [THJ HH] it was suggested to use q 2 as one of the variables. 
The shift of the position of the nucleon pole m m — m would be the unknown 
parameter to be determined from the SR equations. On the other hand 



with U being the single-particle energy of the nucleon. This is the very char- 
acteristics which enters the equation of state. To separate this singularity 
from the other ones a proper choice of the second variable is needed. Con- 
sidering the nuclear matter as the system of A nucleons with momenta Pi, 
introduce 



with p = in the rest frame of the matter. Under the choice of s — 
(p + q) 2 =const we avoid the singularities connected with the excitation 
of two nucleons |IE!~[2IIj- The constant value of s should be fixed by the 




(20) 
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condition that the probe proton is put on the Fermi surface of the matter. 
In the simplified case when the Fermi motion of the nucleons of the matter 
is neglected, we can just assume s = 4m 2 . 

Thus we shall write the dispersion relations for the function G m (q 2 , s) — 
G(q 2 ) with the functions G and G m defined by Eqs. (2) and (20). It was 
shown in j2I] that the nucleon pole is still the lowest lying singularity of the 
function G m (q 2 , s) until we do not include the three- nucleon interactions. All 
the other singularities are lying at larger values of q 2 being quenched by the 
Borel transform due to Eq. (17). Thus we use the "pole + continuum" model 
for the spectrum of the function G m (q 2 , s). 

The OPE coefficients of G m (q 2 ,s) are the in-medium expectation values 
of QCD operators. Thus to use the SR equations one should find the density 
dependence of these condensates. It is not clear "a priori" if the OPE series 
converges indeed. 

3.2 Lowest order OPE terms 

In the lowest order of OPE the expectation values of the lowest dimension 
are involved only. In particular, there are the scalar expectation values 



with qi standing for "u" or "d" quark. There are also the vector expectation 
values 



densates k have nonzero values in vacuum, while the vacuum values of the 
vector condensates vanish, i.e. t>i(0) = 0. Due to the conservation of the 
vector current we find immediately 

(p) , (") 

v (p) = p (25) 

with nf standing for the number of the valence quarks of the flavour "i" 
in the proton (neutron). 

The correlator G m (q 2 ) contains three structures being proportional to 
IfiQ^i an d I with / standing for the unit 4x4 matrix. Thus we 
obtain three QCD sum rules. There are four independent unknowns to be 



K\p) 



(M\q iqi \M) 



(23) 



taking the form v l (p) 
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determined from SR equations. These are the three parameters of the nucleon 
which are the vector and scalar self-energies S s and E s and also the shift of 
the value of the residue A 2 ^ — A 2 . The shift of the position of the pole can 
be expressed through the self-energies, i.e. m m — m = T, v + E s . One more 
unknown parameter is the shift of the position of the threshold — W 2 . 

The explicit form of the SR equations is presented, e.g. in [21]. The 
matching of the rhs and lhs can be achieved in the same interval of the values 
of M 2 as in the case of vacuum. The shift of the position of the nucleon pole 
was found to be a superposition of the vector and scalar condensates 



with the last equality coming form Eq. (21), k = k u + K d . This provides the 
simple picture of formation of the value of m m . Our probe proton exchanges 
quarks with the sea of qq pairs which differs from that in vacuum. This forms 
the Dirac effective mass m* = m + S s . This mechanism is described by the 
first term in rhs of Eq. (26). The exchange with the valence quarks adds the 
second term. 

In the next to leading order of OPE the gluon condensate g(p) = (M\ ^-G^ U \M) 
should be taken into account. We shall see that numerically it is not very 
important. 

3.3 Gas approximation. The role of irN sigma-term in 
nuclear physics 

We saw the lowest OPE to contain the condensates v (p), n(p) and g(p)- While 
the vector condensate is exactly linear in p, the condensates n(p) and g(p) are 
more complicated functions of density. We start with the gas approximation 
in which the matter is treated as ideal Fermi gas of the nucleons. Thus our 
probe proton interacts with the system of non-interacting nucleons. In this 
approximation ^H] 



The matrix elements in rhs of Eqs. (27) and (28) can be related to the 



m m = CtK(p) + C 2 v(p) 
m m -m = U{p) 



(26) 



K(p) 
9(p) 



K(0)+p(N\qq\N) , 
g^)+p{N\^-G^G ia/ \N) . 



(27) 
(28) 
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observables. In particular, 

(N\qq\N) = 2 ° (29) 
m u + m d 

with cr being the pion-nucleon sigma term which is connected to the pion- 
nucleon elastic scattering amplitude T by the relation T = —cr/f% [221 • How- 
ever in the latter relation T denotes the amplitude in certain unphysical point. 
The experiments provide the data on the physical amplitude T ph = —T>/f% 
with S = (60 ±7) MeV. The method of extrapolation of observable amplitude 
to the unphysical point was developed by Gasser et al. [2B]. They found 



a = (45 ± 7) MeV . (30) 

Note that from the point of chiral expansion, the difference £ — a is of higher 
order, i.e. (S — cr)/a ~ m^. 

Note also the physical meaning of the expectation value of the operator 
qq averaged over a hadron state. Anselmino and Forte |2*4^ l2*H] showed that 
under reasonable model assumptions it can be treated as the total number 
of quarks and antiquarks. 

As to gluon condensate, the expectation value is 



(N^G^G^N) = -~(m-£m,-(A%g,|A0) (31) 

with j standing for u, d and s quarks. This equation comes from the averaging 
of the QCD Hamiltonian over the nucleon state with the account of the 
additional relations found in (2H1 . In the chiral limit only the strange quarks 
contribute. In the chiral SU(3) limit the second term in brackets in rhs of 
Eq. (31) turns to zero. 

Solving the SR equations in the leading order of OPE in the gas approx- 
imation we find the potential energy to be [201 



U{p) = bigg[mv{p) - 32(«(p) - k(0)) 



GeV' 2 (32) 



with the difference n(p) — k(0) being described by the second term of the rhs 
of Eq. (27). Thus the potential energy is presented as the superposition of 
a positive term proportional to the vector condensate and a negative term 
proportional to the scalar condensate. At the saturation point p = p = 
0.17 fm -3 = 1.3 ■ 10 _3 GeV 3 we find the two terms in the rhs of Eq. (32) 
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to be 200 MeV and -330 MeV. The gluon condensate adds about 10 MeV 
to the vector term in the next to leading order of OPE Similar results were 
obtained in |2Zj in another SR approach based on the dispersion relations in 
go- 

These are the common points between the SR approach and the Walecka 
model [213 ■ Note, however, that we did not need the fitting parameters like 
g u and g a of QHD. The exchanges by the strongly correlated quarks (i.e. 
by the mesons) are expressed through the exchanges by the uncorrelated 
quarks. The interactions in the vector channel are calculated explicitly. The 
interactions in the scalar channel are expressed through the observable ttN 
sigma-term. Hence, in SR approach the cr-term determines the linear part of 
the scalar interactions. 

3.4 Beyond the gas approximation. A possible satura- 
tion mechanism 

Now we shall try to go beyond the gas approximation, remaining, however, in 
the lowest orders of OPE. Account of the interactions between the nucleons of 
the matter does not change Eq. (25) for the vector condensate. However the 
scalar condensate obtains additional contributions S(p) caused by averaging 
of the operator qq over the meson cloud. Thus we have 

k( P ) = K {0) + p(N\qq\N)+3(p) (33) 

with a nonlinear behavior of S(p) at small values of p. 

Assume that the meson cloud consists of all kinds of the mesons (tc,uj, 
etc.). It was shown in JH] that in the chiral limit m 2 — > (neglecting also 
the finite size of the nucleons) one can obtain the function S(p) as a power 
series in Fermi momenta pp ~ p 1//3 . The lowest order term ~ ppp comes 
from the one-pion Fock term (known also as the Pauli blocking term). The 
two-pion exchanges with the nucleons in the intermediate states contribute 
as pp F |28j . The heavier mesons contribute as ppp ~ p 2 . 

Even beyond the chiral limit we expect the pion cloud to provide the lead- 
ing contribution to the nonlinear term S(p). This is because the contributions 
of various mesons X contains the meson expectation values (X\qq\X) = nx 
with nx being the total number of quarks and antiquarks. We can expect 
nx ~ 2. However the pion expectation value calculated by the current al- 
gebra technique provides (7i\qq\n) = — m + md ~ 12. Thus the pion cloud 
contribution is enhanced. 
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Note that the simplest account of the nonlinear terms in the scalar con- 
densate signals the possible saturation mechanism. Assuming the chiral 
limit vr? v = and including the Pauli blocking term only we can present 
k(p) = «(0) + ;^^P - 3-2^p EH with p F0 = 268 MeV/c being the Fermi 
momentum corresponding to the saturation value of density po- This provides 
the potential energy 



U(p) 



198 - 42 

m u + m d 




MeV , (34) 



which contains the S-term and the pion-nucleon coupling constant g^NN as 
the only parameters. After adding the kinetic energy the energy functional 
obtains the minimum at p = p if we put £ = 62.8 MeV (a = 47.8 MeV) 
which is consistent with the experimental data. The binding energy appears 
to be e = —9 MeV. The incompressibility coefficient which defines the shape 
of the saturation curve also has a reasonable value K m 180 MeV. 

Of course, the results for the saturation should not be taken too seriously. 
They are very sensitive to the exact value of the £ term. This is caused by 
the simplified model of the nonlinear effects. The rigorous treatment of the 
pion cloud requires the account of the multinucleon effects in the propagation 
of the pions [301 EI]- Inclusion of these effects jHSISS] still provides S(p) < 
at the densities close to the saturation value. Thus the nonlinear behaviour 
of the scalar condensate may be responsible for the saturation properties of 
the matter. However the results of this subsection can be considered only as 
the sign that further development of the approach may appear to be fruitful. 

Thus the nonlinear behaviour of the scalar condensate is a possible source 
of the saturation. Here we find a certain analog of the QHD saturation mech- 
anism. Recall that in Walecka model the saturation caused by a complicated 
dependence of the "scalar density" of the nucleons ps = J (f^i^^G ^ ~ P) 

on the density p = J j0sO(p F - p). 

Anyway, further development of the approach requires investigation of the 
higher order OPE terms as well as the analysis of the condensates beyond 
the gas approximation. 



3.5 Higher order OPE terms 

As we have seen, the vector and scalar condensates only contribute in the 
leading order of OPE. The gluon condensate contributes to the terms of the 
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relative order q~ 2 of the OPE of G m (q 2 , s). We saw it to be numerically small, 
as well as some other contributions of this order which we did not discuss 
(see [33J). However there are some problems with the terms of the order q~ 4 
which contain the four-quark condensates h XY (p) = (M\qT x qqT Y q\M) . 

The expectation values of the four-quark operators are not well estab- 
lished even in vacuum. The usual assumption is the factorization approxima- 
tion |H] with the intermediate vacuum states dominating in all the channels. 
As to the in- medium values of the scalar condensate, one can separate the 
configuration with one of the products qq acting on vacuum states while the 
other one acts on the nucleon states [20] • Thus in the gas approximation 

h ss (p)-h ss (0) = 2(0\qq\0)(N\qq\N)p+ (N\(qqqqU\N)p . (35) 

In the second term of the rhs all the quark operators act inside the nucleon. 
Equations (6) and (29) enable us to find the magnitude of the first term of the 
rhs. If the first term estimates the values of all the four-quark condensates 
indeed, the term of the order q~ A of OPE is numerically larger than the 
leading term at the values of the Borel mass where the SR equations are 
solved. This would cause doubts in the convergence of OPE. Fortunately the 
situation is not as bad as that. Celenza et al. |33] demonstrated that there 
is a strong cancellation between the two terms in rhs of Eq. (34). However 
the lack of information on the four-quark condensate have been an obstacle 
for the further development of the approach during many years. The recent 
calculations |33| are expected to improve the situation. 

3.6 New knowledge 

Independently of the magnitude of the contribution of the four-quark con- 
densates, the SR predict certain new features of the nuclear forces. 

3.6.1 Anomalous structure of the nucleon-meson vertices 

As we have stated earlier, the QCD sum rules can be viewed as a con- 
nection between exchange of uncorrelated qq pairs between our probe nu- 
cleon and the matter and the exchange by strongly correlated pairs with 
the same quantum numbers (the mesons). In the conventional QHD picture 
this means that in the Dirac equation for the nucleon in the nuclear matter 
(q — V)ip = (m + $)?/> the vector interaction V corresponds to exchange by 
the vector mesons with the matter while the scalar interaction $ is caused 
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by the scalar mesons exchange. In the mean field approximation the vector 
interaction V is proportional to density p, while the scalar interaction is pro- 
portional to the "scalar density" p s which is a more complicated function of 
the density p. Thus V = V(p), while $ = $(p s ). We have seen that QCD 
sum rules provide similar picture in the lowest orders of OPE: vector and 
scalar parts G V;S of the correlator G m depend on vector and scalar conden- 
sates correspondingly: G™ = G" l (v(p)); G™ = G™(k,(p)). However we find 
a somewhat more complicated dependence in the higher order OPE terms, 
say, G™ = G™(k(p),v(p)), depending on both scalar and vector condensates. 
This is due to the four-vector condensates. In particular, the scalar-vector 
condensate {M\uud n fod\M) contains the contribution (0\uu\0) (M\djod\M) . 
This term is proportional to the vector condensate, contributing, however, 
to the SR for the scalar structure. In QHD picture this corresponds to the 
explicit dependence of the scalar interaction on the density p. Such de- 
pendence is not included in QHD, at least on the mean field level. These 
contributions correspond to the anomalous structure of vertex of the inter- 
action between the nucleon and the scalar meson in the meson exchange 
picture. Similar situation takes place for the condensate (0\uu\0) (M\uu\M) 
contributing to the SR for the vector structure. 

3.6.2 Charge-symmetry breaking forces 

There is an old problem of the difference between the strong interactions of 
the proton and neutron with the systems containing equal numbers of pro- 
tons and neutrons [36j. In framework of SR method the difference can be 
caused by the explicit dependence of the function G m on the current quark 
masses and on the isospin breaking expectation value 7 = ^m\u,u\m)^ • ^ ne 
problem was attacked by the SR approach in a number of works |HZj. The 
proton-neutron binding energy difference was expressed through the quark 
mass difference — m u and the condensate 7 under various additional as- 
sumptions. The SR analysis provided some qualitative results which may be 
useful in the building of the charge-symmetry breaking nuclear forces (CSB). 
One of them is the importance of CSB in the scalar channel. Earlier there was 
a common belief that the CSB in the vector channel are responsible for the 
effect. Another point is that the mixed structures described in Subsec. 3.6.1, 
manifest themselves in the leading terms of OPE. Thus explicit dependence 
of vector and scalar forces on both p and p s may become important. 
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3.7 A self-consistent scenario [32,33] 

A rigorous calculation of the contribution of the pion cloud to the scalar 
condensate n(p) requires the consistent treatment of the pion dynamics in 
nuclear matter. The nonlinear term S(p) can be presented as 

s( p ) = 2m " ^(p) (36) 

rn u + m d dm\ 

with being the pion contribution to the nucleon self-energy. This can be 
obtained by using the presentation S = dV/dm q (with V standing for the 
interaction energy) obtained by Cohen et al. [HE] and taking into account 
the pion contribution only. The self-energy contains the in-medium pion 
propagator D which satisfies the Dyson equation 

D = A) + D UD • (37) 

The propagator D differs from the free propagator D due to the particle- 
hole excitations described by the polarization operator II. The latter can be 
expressed through the amplitude of the forward ttN scattering (this provides 
the value n ). The short-range correlations can be described by means of the 
Finite Fermi System Theory (FFST) introduced by Migdal [SHI- If only the 
nucleon-hole excitations are included, the operator n turns to II = Il /(1 + 
PAfllo) with gw being a FFST constant determined from the experimental 
data. The account of long-range correlations [35] modifies the effective value 
of g N . 

The function S(p) thus depends on the effective mass of the nucleon m*(p) 
and on the ttN coupling constant g* NN (p)- The latter can be presented 
through the fundamental parameters /„- and axial coupling constant g& by 
Goldberger-Treiman relation jJU] 

9nNN 9A / qc x 

"2^T ~ 2U ' (38) 
which can be expanded to the case of the finite density. 

On the other hand, the SR provide the dependence of the nucleon effective 
mass on the condensate n(p). If we succeed in describing the dependence of 
9a(p,k(p)) (the first steps were made in [20]) and of f n (K(p),p), we shall 
come to the set of self-consistent equations 

* %) 
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The first three equations should originate from the sum rules. The last one 
is just the combination of Eqs. (33), (36) and (38). One should add similar 
equations for the delta-isobar. 

This would enable to describe the baryon parameters in the hadron phase 
of the baryon matter up to the point of the chiral phase transition where 

K{p) = 0. " 

3.8 A sub-plot: Goldstone pions never condense [32,33,41 

The possibility of the "pion condensation" was first discussed by Migdal 
[7j. The observation is that at certain value of density p = p v the pion 
propagator presented by Eq. (36) obtains the pole at the energy e w = 0. 
This would signal the degeneracy of the ground state of the system. The 
ground state contains an admixture of the oscillations with the quantum 
numbers of the pions. The value of p v appeared to be sensitive to the values 
of FFST constants, being p v > 2p in most of the hadronic models. 

The analysis shows that the singularity of the pion propagator in the 
point of the "pion condensation" leads to the divergence of the function S(p) 
expressed by Eq. (35). This provides n(p) — > +oo at p — > p n , while the 
density increases. On the other hand, k(0) < 0. Thus k = at certain 
point p c h between the zero value and p n , i.e. < p c h < p w . However, 
since n(p ch ) = 0, the chiral symmetry is restored at p = p c h- One cannot 
expand the hadronic physics of the densities close to the saturation value 
p to the region p > p c h. Thus the "pion condensation" point cannot be 
reached in framework of the existing hadronic models. Anyway, once the 
chiral symmetry is restored, the pions do not exist as the Goldstone bosons 
any more. Thus, there is no "pion condensation" of the Goldstone pions. 

The subject is analyzed in details in the papers |^TJ|32]- Although not be- 
ing connected with the SR directly, this result is the outcome of investigation 
of the scalar condensate stimulated by studies of the sum rules. 

4 Summary 

We saw that the first steps in the application of QCD sum rules method to in- 
vestigation of the in-medium nucleon parameters are successful. The effective 
mass and the single-particle potential energy where expressed through QCD 
condensates in the lowest OPE orders. The quark and gluon condensates 
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are expressed through the observables. The pion-nucleon a-term appeared 
to determine the value of the scalar forces. The approach reproduced one 
of the key points of QHD picture: the nucleon moves in superposition of 
scalar and vector fields which cancel to large extent. The numerical values 
are consistent with those of QHD. However the approach does not use the fit- 
ting parameters and avoids the controversial conception of the heavy meson 
exchange by the point nucleons. 

The approach provides also some new knowledge about the nuclear forces. 
These are the anomalous structures of the meson-nucleon vertices. Such con- 
tributions usually are not included in QHD. Another point is the importance 
of the scalar channel in the charge-symmetry breaking forces. 

Even the simplified model for the in-medium scalar condensate provides a 
possible mechanism of saturation. It is caused by the nonlinear contribution 
of the pion cloud to this expectation value. In the rigorous treatment of the 
pion dynamics the sum rule for the effective mass of the nucleon m* and 
the expression for the pion contribution to the condensate form the set of 
self-consistent equations. The complete set of the equations should include 
the sum rules for the axial coupling constant g* A in medium and for the in- 
medium pion decay constant /*. Investigation of the complete set of the 
equations is the subject of the future work. 

The development of the method requires calculation of the troublesome 
four-quark condensates. The calculation of complete set of these expectation 
values in framework of the convincing models is in progress. 
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